Supersymmetry and fluctuation relations for currents in closed networks 
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We demonstrate supersymmetry in the counting statistics of stochastic particle currents and use 
it to derive exact nonperturbative relations for the statistics of currents induced by arbitrarily fast 
time-dependent protocols. 



Introduction. The discovery of fluctuation theorems 
and nonequilibrium work relations pQ has stimulated con- 
siderable interest in nonequilibrium statistical mechanics 
and theory of counting statistics [2]. Fluctuation the- 
orems apply to evolution of thermodynamically impor- 
tant characteristics, such as work, dissipated heat and 
entropy. The utility and even a proper definition of these 
thermodynamic concepts in the framework of mesoscopic 
nonequilibrium physics are still a subject of considerable 
debates. Hence, it is important to obtain exact relations 
that do not directly rely on the thermodynamic con- 
cepts, but rather describe unambiguous characteristics, 
such as statistics of particle currents in systems driven 
by time-dependent fields. Progress in this direction has 
been relatively modest. Fluctuation theorems have led 
to exact relations for statistics of particle currents only 
in nonequilibrium steady states |2J [3] , or for driving pro- 
tocols that do not break time-reversal symmetry [4j. 

In this letter, we present exact relations for statistics 
of currents in strongly driven mesoscopic stochastic sys- 
tems. Being akin to known fluctuation theorems, a part 
of our exact result is not directly related to the condi- 
tion of microscopic reversibility but rather follows from 
supersymmetry of the counting statistics of currents. 

Model. In this section, we introduce a notation that 
will be used to formulate our results. Consider a graph 
X = (Xq,Xi), with X and X\ being the node and link 
sets, respectively, and a Markov process describing a par- 
ticle motion on X with a set of transition rates, kfj, 
of jumping from the node i to the node j through the 
link a. Each link connects two different nodes. Fig. [T] 
shows an example of such a graph with 6 links and 5 
nodes. To define positive direction of currents, we pre- 
scribe arbitrary orientations on all links. We denote 
by Co(X) and C\{X) the vector spaces of distributions 
P = (pj\j ^ -^o) and currents on links J = (J a \a 6 X%), 
respectively. Nodes and links are labeled by Latin and 
Greek indices, respectively. We also introduce a notation, 
da = = (doa, diet), which means that da is the ori- 
ented border of the link a, and it can be represented by 
an ordered pair of nodes where j is a node to which 
the link a points and i is the node from which this link 
originates. With a minimal abuse we will use the same 
notation, da = {i,j}, for non-oriented boundary. It is 
convenient to view a current as a set of components (J z "), 
with {i,j} = da, such that = — J™. Conservation of 




FIG. 1. A closed graph with 5 nodes and 6 links, repre- 
senting a Markov chain with five states. Links represent al- 
lowed transitions with rates kfj = ke^^- Ei ~ Wa+ii i\ Nodes are 
characterized by energies of wells Ei, i = 1,...,5. Links 
are characterized by sizes of barriers W a , by parameters 
F a = F% = ln(k™/k«) = f-j- ffi + Ei - Ej, and by count- 
ing parameters \a = X?j that are used to count number of 
transitions through the link a. Here a — 1, . . . , 6. 

particles requires that J^ - pj = 1. The evolution of the 
probability vector is given by the master equation, 

dtp = Hp, (1) 

where we will call H the master operator. Let \j) and (i\ 
be the bra- and ket- vectors over the space Cq{X) with the 
only nonzero unit components at j-th and i-th positions 
respectively. In this basis set, Hij — Y^*~ ^'"^ f° r 
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Detailed balance (DB) is not assumed, the deviation 
from DB is quantified by the entropy function (EF) that 
associates F(l) = ^Y\ ae i(kfj / k^) with a closed path 
I, so that, for time-independent parameters, e F ^ repre- 
sents the ratio of the probability of a stochastic trajectory 
I to its time- reversed counterpart. Here we introduced 
a notation, Y\ a( zi{kfj / k^) , which means the product of 
kfj/kj 1 } over all links of the loop, where we assume an 
arbitrary direction of motion along the loop as positive 
and take = da if the link a points along the positive 
loop direction and (j, i) = da if the link points against 
direction of the loop. Later, we will use analogous nota- 
tion for the sum over loop links, Ylael' Note that DB 
means that F(l) = for any closed path I. We will con- 
sider only driving protocols that conserve the entropy 
functions F(l) at any cycle I of a graph. 
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We focus on the case of periodic driving, when the ki- 
netic rates depend on time in a periodic manner, accord- 
ing to a driving protocol with the time period r. Such 
a steadily driven system eventually enters a regime with 
periodically changing population probability vector, i.e. 
p(t + t) = p(t). We can introduce the currents per pe- 
riod of the driving protocols J c = (f™ T dtJ{t))/n. Since, 
in the limit n — > oo, particles cannot accumulate in any 
node, the sum of currents entering any node is equal to 
the sum of the currents leaving this node. We will call 
currents with this property the conserved currents. 

The entropy function is naturally extended to a linear 
entropy functional on currents, given by 

£ (J c ) = ]T F a J c a , (2) 

where we introduced the vector F, whose components 
are indexed by the links, given by F a = F°j — — F^ = 
\n(k°j / k^) . It is important to note that even though 
the parameters kfj, and hence F£*, are time-dependent, 
the entropy functional £ defined on conserved currents 
J c is time- independent, provided all F (I) are time- 
independent. One can see this, e.g., by noticing that 
there is a basis set in the space of possible conserved 
currents, which consists of constant unit-valued currents 
circulating in each independent cycle of a graph and hav- 
ing zero values on all other links. The contribution of 
each such independent conserved current to the entropy 
functional is just the entropy function for a correspond- 
ing loop of a graph, and any conserved current is just a 
linear combination of such circulating basis currents. 

Fluctuation relation for currents (FRC). Consider a 
Markovian kinetics of a particle on a graph, with constant 
entropy function, F(l), and periodically time-dependent 
kinetic rates, kfj. The probability distribution of con- 
served currents, generated per period of driving, has 
the large-deviation (LD) form P(J c ,n) ~ e nS{J a ) ^ w ^ 
S(J C ) being referred to as the LD function. Consider the 
following symmetry property of the LD function S(J C ) 
of the conserved currents 

S f (J c ) -S b (-J c ) = £{J C ), (3) 

where / and b stand for the original (forward) and time- 
reversed (backward) driving protocols respectively. Ki- 
netic rates in backward and forward protocols are related 
by k% b (t) = kf/{r - t). Eq. (B is known to hold if the 
kinetic rates are independent oftime 3, 4 J but it does not 
hold for general time-dependent rates [I]. Nevertheless 
we can show, and this is the content of our FRC, that 
Eq. ([3| does hold for two types of periodic driving pro- 
tocols, for which, in addition to conservation of entropy 
function, either of the two conditions is satisfied: 

(i) the ratios kfj/k^ of the forward/backward rates 
are kept time- independent at all links a. This can be de- 
scribed by introducing time dependent parameters, W ai 
on the links, such that kfAt) = k%e w "^ . 
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(ii) the branching ratios kfj /kfj, are kept time- 
independent at all nodes. This can be described by in- 
troducing time-dependent parameters, Ei, on nodes, such 
that fcg(t) = fcge- JS( M. 

It is possible to write kinetic rates in the form kfj = 

e Wa ~ E * + f*i , with a constant vector / = (f°j\a G X\) on 
the links, and, generally, / 2 " ^ /?• . For the case of a 
network that describes transitions of a physical system 
among deep free energy minima in its phase space, pa- 
rameters W a , Ej, and ff. have a clear physical interpreta- 
tion [7]. They correspond, respectively, to the potential 
barrier separating metastable states along the path a, 
to the size of the energy of a well in the node i, and to 
the effect of an external force acting on the system along 
the path a. Adopting this terminology, we can formu- 
late the FRC as the statement of the validity of Eq. ^ 
for time-dependent protocols in which either (i) only po- 
tential barriers are driven or (ii) only node energies are 
driven. 

It is important to note that parameterization of ki- 
netic rates by the set (E, W, f ) is not unique since it 
is possible to redefine parameters f°j — > — w a — Si, 
W a — > W a + w a , and Ei — > Ei — Ei, and obtain the same 
set of kinetic rates kfj. This change of parameters pre- 
serves the entropy function. In fact, a set of kinetic rates 
is fully determined by the entropy function F(l) and the 
sets E = (Ej\j G X ) and W = (W a \a G Xi). The 
entropy functions for any loop I can then be expressed 
in terms of / alone, i.e. F(l) = J2 a ei(f? 3 - fji); this 
means that the entropy function and the entropy func- 
tional ([2]) are invariant with respect to the rate trans- 
formations kfj 1 — y kfje Wa ~ £i . We will use this property 
in our derivation of FRC, namely, if we can prove FRC 
for some choice of the vector / which corresponds to a 
given form of the entropy functional, then FRC is valid 
for any other choice of a constant / that corresponds to 
the same entropy functional. During the derivations of 
(i) and (ii), different choices of / will be used to apply 
the symmetries of the problem. 

Operator derivation of case (i) and Lagrangian inter- 
pretation. We will derive FRC by considering the sym- 
metries of the Legendre transform, w x , of S(J C ), where 
X is the variable conjugated to J c . In the literature, 
lu x is often called the cumulant generating function of 
currents [5]. An operator approach to deriving the cu- 
mulant generating function of conserved currents is based 
on introducing the twisted master operator, H x , parame- 
terized by the multi- variable argument x = (Xal a € X\) 
(or alternatively by a set of antisymmetric components 
Xfj = ~Xfn called counting parameters) of the gener- 
ating function. Here the word "twisted" reflects the 
way the operator H x is constructed, i.e. by multi- 
plying (twisting) off-diagonal elements of H by corre- 
sponding e Xi i factors, namely, the operator H x is ob- 
tained by replacing the off-diagonal components of H 
with (i\H x \j) — e%ji ^i! an( A with the diagonal 

components remaining the same [5]. 
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FIG. 2. Equivalence relations among counting parameters. 

After many periods of driving, lo x is determined by the 
largest eigenvalue, e Wx , of the evolution operator 



U x = Texp U H x (t)\ , U x \i/> x ) = e^ x ). (4) 

Eq. Q is then equivalent to the following symmetry 
property of lo x : 

u) L= u -x-f ( 5 ) 

For case (i), Eq. ^ can be derived in a concise way. 
By direct verification, we have the symmetry 



x-f, 



(6) 



where H x is the transpose of H x and, by the condition 
imposed on the driving protocol, (i), the vector F is time- 
independent. The transposition changes the ordering of 
operator products, while maintaining the eigenvalues un- 
changed. This implies {U{) T = U b _ x _ F . Ec t- ® follows 
from the fact that for a constant set E we can always 
redefine the set / so that f = F. 

Our elementary derivation of (i) is complemented by 
a simple interpretation in terms of probabilities of parti- 
cle trajectories, P(x), where a closed particle trajectory 
x — (I, t) is represented by a closed path I on X, and tem- 
poral data t (the times when the jumps occurred). We 
then have the symmetry property pf(x)/P b (x r ) — e F ^ 
for the probabilities of the original trajectory and time- 
reversed trajectory x r for the time-reversed protocol. 
For case (ii), the above arguments cannot be applied, 
since F depends on energies, which in case (ii) are time- 
dependent, so that the symmetry |6]) cannot be consid- 
ered equivalent to ([5| . Neither does case (ii) have a sim- 
ple interpretation in terms of stochastic trajectories. We 
will derive Eq. ([3| for case (ii) using an additional super- 
symmetry property of the evolution with twisted master 
operator. 

Super symmetry for master equation. A hidden super- 
symmetry of Langevin dynamics and motion on cyclic 
graphs has been discussed in the literature (see e.g. 
in [5]). Here we show that a hidden supersymmetry 
can be found in Markovian evolution on an arbitrary 
graph. More importantly, this supersymmetry can be 
extended to make it a property of the counting statis- 
tics of currents. The master operator has a representa- 
tion H — — QJ, with J being the current operator, and 



a). The operator 
i) for da = (i,j). 



Q acts from C\ to C by Q\a) = \j) 
Here \a) is the vector in C\ with a unit entry correspond- 
ing to link a and zero otherwise. The operator Q plays 
the role of a discrete counterpart of the div-operator [7] . 
Therefore, for conserved currents, we have QJ C = 0. 

We further introduce an operator H = —JQ that acts 
in the space C± of currents, and notice that obviously 
HQ = QH. Viewing H = H © H and Q as operators 
acting in C = C © Ci, we have [Q, H] —Q. Considering 
Co and C\ as even and odd components, of the (super) 
space, Q becomes an odd operator that commutes with 
the master (super) operator H, and H can be written as 
the anticommutator of — Q and J, which closes the al- 
gebra; therefore, the term supersymmetry is totally ap- 
propriate here. Supersymmetry connects evolution of the 
probability distribution and current distribution compo- 
nents of the superspace, C, thus allowing the standard 
master equation to be reformulated in terms of its coun- 
terpart that describes evolution in the space of currents. 
Our key observation for proving (ii) is that supersymme- 
try for master operator also holds for the twisted oper- 
ator, H x . By reformulating the problem of finding the 
generating function in terms of the superpartner of H x , 
we will show that a derivation of case (ii) becomes as 
simple as for case (i). 

Supersymmetry for twisted operators. To identify su- 
persymmetry on the level of twisted operators we need to 
come up with a procedure for twisting the current mas- 
ter operator H. Since the off-diagonal elements of H are 
between the links that share a common node, it is rea- 
sonable to represent the twisting data by C = (CcUli € 



da D dj) with antisymmetric components C, 3 a 
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This allows a family of twisted current master oper- 
ators to be introduced by 



9j=(j,m) 9j=(m,j) 
+ ( J2 e ^ k ]m- E e< ^ k J™ 

d^f— (m,z) d r y—(i,m) 

+ E eCfc 7 m - E eC - fc 7Jl7>, (7) 

for a = as well as a family of twisted evolution 

operators, using a definition, similar to Eq. Q. 

Parameters, x an d C have a purpose of counting how 
many times a particle passes, respectively, through links 
and nodes of a graph along specified directions. The 
number of independent conserved currents, however, is 
smaller than the sizes of these vectors. The informa- 
tion about conserved currents is contained in depen- 
dence of the generating function only on expressions, 
r x (0 = E„ 6 iX§ or T c (l) = £ jei a 7 , where the last 
sum runs over nodes that belong to the loop i, and in- 
dexes a and 7 in correspond to the links that, respec- 
tively, precede and follow the node j along the positive 
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loop direction. We will call two sets, x an d C equiva- 
lent, if r x (i) = r^(Z) for any closed path, I. A simple, 
yet important observation is that x 1S equivalent to C if 
and only if there is a set cp = ((fj a \j £ da) such that 
X?j = ¥>ja - Via and £q 7 = Vja - <Pj~/- We additionally 
illustrate the meaning and relations among the counting 
parameters in Fig. [2} It is straightforward to verify that, 
provided <p establishes equivalence between x an d Cj the 
following supersymmetry relation takes place: 

H X Q V = Q V H C , U X Q V = Q V U C , (8) 

where Q v is the twisted supersymmetry operator, de- 
fined by Qip\a) — e v i"\j) — e Via \i) for a link a = 
Although Eq. ^ can be verified directly by 
inspecting the matrix elements, it is instructive to 
note that it follows immediately from the represen- 
tations H x = —QcpJcp, and = —J V Q, V , where 
J v is the twisted current operator, given by Jtp\i) — 

(J2 d a 7 m fc 5 e ~^ " T^= {1 ' J) Age"*"-) \a). We note 
also that Q v is time-independent so that the evolution 
operators with U x and have the same sets of nonunit 
eigenvalues even if the parameters E are periodically 
driven since if \ip) is the eigenstate of the U^, then Q v \ip) 
is the eigenstate of U x . 

Derivation of case (ii). To reboot the operator deriva- 
tion of Eq. ([5]) for case (ii) , we introduce the superpart- 
ner, F — (F^\j £ daHdj), of F with the antisymmetric 
components = -F^ a = — ln(fc^/A:J m ). Obviously, 

the entropy functions satisfy F(l) — F(l) = J2jei ^7 
for any closed path I, referred to as consistency. This 
means that F represents a vector which can be chosen to 
be equivalent to /. One verifies directly that the operator 
|7| satisfies the relation, = H^g, and note that in 

case (ii) the data F is time-independent, which results 



in (U£) T = U b ^ ~ and further in the symmetry relation 
ujj* — uj b ~ for the eigenvalues of the evolution opera- 
tors in the current space. We further naturally choose C 
equivalent to xi this results in — £ — F being equivalent 
to —x ~ /; due to consistency between / and F. This 
allows us to apply supersymmetry [Eq. Q], which results 
in uj x — uj^ and u h _ x _^ = uj b ^ ~, and further in Eq. (JiJ), 
which completes our derivation of FRC. 

Discussion. The FRC is applicable to the case when 
either energies of discrete states or the heights of barriers 
that separate the states are varied. This regime can be 
realized for catenane molecular motors or electric circuits 
in incoherent regime, when the gate voltages of the quan- 
tum dots are varied in time [7] . Statistics of currents can 
be probed in single molecule measurements [S] or in ex- 
periments with nanoscale electric circuits [5]. The FRC 
demonstrates that there are fluctuation theorems that do 
not directly follow from the relations between the prob- 
abilities of forward and time-reversed trajectories. We 
expect that the FRC is only one of many possible appli- 
cations of the supersymmetry of the counting statistics 
of currents. 
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